A cyclic cooling algorithm 
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We introduce a scheme to perform the cooling algorithm, first presented by Boykin et al. in 2002, 
for an arbitrary number of times on the same set of qbits. We achieve this goal by adding an 
additional SWAP-gate and a bath contact to the algorithm. This way one qbit may repeatedly be 
cooled without adding additional qbits to the system. By using a product Liouville space to model 
the bath contact we calculate the density matrix of the system after a given number of applications 
of the algorithm. 

PACS numbers: 03.65.-w, 03.67.-a, 03.67.Lx, 05.30.-d 



Algorithmic cooling is a method to obtain highly 
polarized spins in a spin system without cooling down 
the environment. It may for example be used for medical 
magnetic resonance imaging to improve the resolution by 
cooling down a subset of nuclear spins of a patient with- 
out cooling of the patient himself, or for the preparation 
of the ground state of a quantum computer by means of 
the computer itself, that means that no external cool- 
ingmechanism would have to be attached to the system 

m 

The spin to be cooled down (a nuclear spin for exam- 
ple) has to couple weakly to the environment. In addition 
one uses some rapid relaxing spins to transport energy 
out of the system. The transportation of energy from 
the cooled spin to the others is achieved in a strictly 
non classical way by applying a quantum algorithm to 
the system, therefore the spins are further referred to as 
qbits. 

Recently, Boykin et al. [|[ have developed a quantum 
algorithm to cool down a single qbit with the aid of two 
auxiliary qbits. Initially the system is prepared in an 
equilibrium state with all spins at the same inverse tem- 
perature (3(0). (Note that we label all quantities belong- 
ing to the n th application of the algorithm by (n), so the 
initial state and the temperature of the bath, introduced 
later on, are labeled (0).) By applying several quantum 
gate operations one spin is cooled down by transferring 
energy to the others. The algorithm itself consists of 
a controlled NOT (CNOT) gate and a controlled swap 
gate (CSWAP) [|. The CSWAP is a 3 qbit gate which 
swaps qbit 1 with qbit 3 if qbit 2 is |0), otherwise it 
does nothing. This leads to an increase of the inverse 
temperature (3(1) of qbit 1 to approximately 3/2 /3(0) 
(cf. [1, @). Having applied the algorithm once, the ini- 
tial state is recovered by two further applications of the 
algorithm. However, by cooling down two other qbits by 
applying the same algorithm as described above to two 
additional sets of 3 qbits allows a second application of 
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FIG. 1: Cyclical cooling algorithm: 3 quantum gates are ap- 
plied, first a SWAP gate then a CNOT gate and finally a 
CSWAP gate. Bath contact at inverse temperature /3(0) is 
symbolized by the boxes on qbit 2 and 3. 



the algorithm to the cooled qbit triple with reduced ini- 
tial inverse temperature [3(1). Thus, one qbit could be 
cooled down to the total of 9/4 /3(0). Possessing an un- 
limited number of qbits the method is, in principle, able 
to reach arbitrary low temperatures for a single qbit, but, 
of course, due to an exponential growth of resources 0. 

As an improvement in the same paper Q the use of 
rapidly thermalising qbits is suggested. That means that 
after the application of the algorithm qbit 2 and 3 are 
coupled to a polarization heat bath with inverse tem- 
perature (3(0) and thus relax back to their initial state 
(infinite bath contact time). These two qbits can now be 
used to cool down a second and third qbit to use those 
for a second cooling step for one of the cooled ones. This 
system is not closed any more and so Shannon's bound 
does not apply. But even so the growth of the number of 
needed qbits is slowed down it is nevertheless exponen- 
tial. 

Instead of using an infinite number of qbits to reach 
an arbitrary low temperature it would be highly desirable 
to have an algorithm which reaches at least some lower 
temperature without using more qbits. As will be shown 
below, this is obtained by starting with a SWAP gate 
between qbit 1 and 3. After the SWAP again a CNOT 
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and CSWAP gate is applied to the system in the same 
manner as in Boykins algorithm. Now one has to wait 
until the auxiliary qbits 2 and 3 are relaxed back to the 
initial temperature by the coupling to a heat bath (see 
Fig. [I]). For this last step we investigate two cases: Total 
relaxation to the bath temperature (infinite contact time) 
and a finite coupling time r to the bath, that means 
coherences in the system are not totally damped and the 
initial occupation probabilities are not entirely regained. 

This new algorithm can be applied for an arbitrary 
number of times to the same set of qbits. We thus call 
this algorithm cyclical, despite the fact that after one ap- 
plication the system does not return to its initial state. 
We are then interested in the final inverse temperature 
/3(n) of the first qbit after the n-th application in depen- 
dence of system parameters and relaxation times between 
subsequent application steps. In principle, the complete 
process could be seen in the context of thermodynamical 
machines (refrigerators) cooling down a finite "environ- 
ment" (here only the single spin 1). 

Considering more than three qbits this cyclic algorithm 
is able to cool down one half of the accessible qbits by 
applying it on sets of 4 qbits, first cooling qbit 1 with the 
aid of qbit 2 and 3, then cooling qbit 4 with the qbits 
2 and 3, again. If one was capable of running quantum 
gates on any combination of qbits, all but two qbits of 
the system may be cooled down. 

The non-interacting 3 spin system is described by the 
Hamiltonian 
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H = Y / ^a z ( f i) (1) 

where AE^ specifies the respective Zeeman splitting of 
qbit p. Each quantum gate is represented by a unitary 
transformation U. Thus we introduce the transforma- 
tion operators [/swap, C^cnot and JTcswap, representing 
the whole algorithm as U = L^cswap^cnot^swap- The 
density matrix after the application of the algorithm, but 
before the next bath contact, is given by 

p { (n) = Ufr(n)U-\ (2) 

where n denotes the number of applications of the algo- 
rithm until now. 

Then the bath contact has to be taken into account 
in order to calculate the final temperature. Since the 
qbits are uncoupled, each one is thermalized separately. 
A standard technique to describe such a bath coupling 
refers to the quantum master equation in Lindblad form 
0, H, H, EH ■ For a single spin p the respective Liouville 
von Neumann equation reads 

+ W Q - t i(2a + p fl a^ - p fl a^a+ - a^a+p^), (3) 

according to the rates W\^q = A/ (l + (l/e)) and Wo^i = 
A/(l + e), e = cxp[A£73(0)] and the bath coupling 



strength A. (Here we used for all numerical investiga- 
tion A = 0.01.) 

The total Liouville von Neumann equation can be rep- 
resented by the superoperator £ acting on operators of 
the Hilbert space H, here the density operator 

h = £/V (4) 

Sorting the entries of an operator O on TL (for example 
the density operator) into a k 2 dimensional vector (with 
k being the Hilbert space dimension), we define "ket" 
and "bra" like vectors O — > \0) and — > (0\ in this 
super space. Their inner product is defined as (A\B) = 
Ti{A^B}, the trace norm of operators in TL. Operators 
D = \A)(B\ actin g on states \0) in the Liouville space 
are defined as (cf. [ll[) 

£)\6) = \A)(B\6) = Tr{&6}A. (5) 

where the superoperator O represents a k 2 x k 2 dimen- 
sional matrix. For a single spin a convenient basis is given 
by the Pauli operators &i with i — {x, y, z, 0} (oo to 
represent identity). Transferring those Hilbert space op- 
erators into the Liouville space as described above each 
superoperator may be expanded as 

O = Y,O ij \a i )fa\. (6) 

ij 

One of the big advantages of this super space formalism is 
the possibility of writing down a superoperator projecting 
an arbitrary state on a solution of Eq. Just like 

in Hilbert space the time evolution operator [the formal 
solution of Eq. is given by 

p(t) = e^pfM = X(t)pt(n) (7) 

with the limit 

lim %{t) = X(oo) (8) 

t — >OG 

defining the complete thcrmalisation superoperator. 
Based on Eq. ([7]) we find for diagonal density operators 
the time evolution operator in terms of the Pauli basis 

S(T)H^>)(<7o|+ e - 2rA \&,){&,\ 

+ {e- 2TX -l) £ j^\\a Q )(a z \. (9) 

This superoperator represents the thermalisation process 
truncated after a time step r. In this case the bath has 
not jet completely thermalised the spin. To extend the 
thcrmalising superoperator to more than one spin we use 
a product Liouville space with the basis \&i) ® |<Tj) <S) 
^k) = \&ijk)- In this basis the respective thcrmalising 
superoperator of qbit 2 and 3 (cf. Fig. [1]) reads 

T 23 (r) = 1®£(t)®X(t). (10) 
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FIG. 2: Trend of /3(n) as a function of cycles n for short 
bath contact time (r = 1/50 Tl). The three curves together 
describe the temperature evolution of qbit 1: (3(n) jumps from 
the upper curve to the lowest, then to the middle one and up 
again. 



For the superoperators of the quantum gates we use 
the corresponding unitary transformation @ on a gen- 
eral density matrix expanded in Pauli matrices, thus, ob- 
taining the superoperator U pU~ l — > ii\p) in the product 
basis. We get the superoperator 05 = £2,3 (t)H for the 
complete algorithm and the state of the system after an 
arbitrary number n of cycles 



\h{n)) = <8 n |/3i(0)). 
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This could be achieved by finding the Jordan decompo- 
sition of 23 and take it to its n-th power. Finally, we get 
the density operator in the Liouville space product basis. 
The result is easily transformed back to the Hilbcrt space 
by using the respective Pauli product basis in the Hilbcrt 
space <Jijk = di ® &j ® &k, by inserting the coefficients 
P{,ijk(n) of \pf(n)) back into the expansion 



ijk 



ijk- 



(12) 



In the case of complete relaxation of the auxiliary 
qbits, i.e., by taking the limit for r — » oo of %2s{t), wc 
have been able to compute the limit n — > oo of Eq. (fTT|) . 
The asymptotic inverse temperature of the cooled qbit 1 
yields 



/3(oc) 



AE 2 + AE 3 
AE 1 



13(0). 



(13) 



It just depends on the Zceman splitting of the qbits in- 
volved in the algorithm and may thus be adjusted pre- 
cisely. 

Unfortunately, for the truncated relaxation, i.e., a fi- 
nite relaxation time r, a complete analytic solution is not 
available, hence we investigate the algorithm numerically 
(with AEi = AE 2 = AE 3 = 1). For short bath contact 
time r the inverse temperature @{n) of the cooled spin 
is shown in Fig. [21 for long bath contact it evolves like 
Fig. Apparently, in Fig. [3] the final stationary state is 



FIG. 3: P(n) as a function of cycles n for long bath con- 
tact time (r = 4T1). Here the final temperature is reached 
approximately after a few steps. 
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FIG. 4: The final inverse temperature /3(300) as a function of 
different bath contact times r. The line at /3(300)//3(0) = 2 
represents the upper limit for infinite bath contact, the line 
at /3(300)//9(0) = 1, 53 marks the inverse temperature of qbit 
1 after one application of the algorithm. 



reached quickly after fewer then ten applications of the 
algorithm. This is not the case for very short relaxation 
times, as can be seen from Fig. [5] 

We measure the time in spin-laticc rcalxation times 
Tl. After each application of the quantum gate we wait 
exactly the same time r to let qbit 2 and 3 relax towards 
the equilibrium temperature of the bath. In Fig. [4] we 
show the final inverse temperature of qbit 1 after n = 300 
applications of the full algorithm (gates plus bath cou- 
pling) in dependence of the bath coupling time r. Note 
that even if we have computed three hundred applications 
of the algorithm for each different waiting time r, the fi- 
nal temperature can already be reached after a couple of 
applications. This is especially the case for t ^> 0.64 Tl 
as could be seen in Fig. [3l 

At r w 0.64 Tl the final temperature is as low as in 
the closed Boykin algorithm. Note, that a single cycle of 
the new algorithm leads to the same temperature as the 
Boykin algorithm. Waiting time r sa 0.64 Tl means that 
qbit 2 and 3 are cooled to 0.77 times the inverse bath 
temperature (3(0) after the first step. 

For t > 0.64 Tl the final temperature is lower than 
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it is in the Boykin algorithm (line in Fig. [4] at w 3/2). 
Thus, the cooling of the algorithm is improved, compared 
to a single cycle. However, a multiple application of the 
algorithm does not always lead to a reduction of the final 
temperature of the cooled spin at least if r < 0.64 IT. 
This can be seen from Fig. [2] as well, where the temper- 
ature after the first application is already smaller than 
after 300 cycles. 

For an easier comparison, we summarize the require- 
ments and final temperatures of the discussed algorithms 
in the following table 



Algorithm 


qbits 


cycles 


final (3{n) 


Boykin 


3 


n = 1 


~ 3/2/3(0) 


Boykin 


3 


n = 2 


= /3(0) 


Boykin unlim. res. 




n = 3 N 


~ (3/2)^(0) 


Boykin+Bath 


2 + 


n = 3^ 


~ (3/2)^(0) 


Cyclic (r =arb.) 


3 


n = 1 


~ 3/2/3(0) 


Cyclic (t = oo ) 


3 


n = oo 


= 2/3(0) 


Cyclic (t = 1/50 Tl) 


3 


n = 300 


= 1.37/3(0) 


Cyclic (t = 4T1) 


3 


n = 300 


= 1.97/3(0) 


Cyclic (t = 4 IT) 


3 


n = 6 


= 1.96/3(0) 



Finally, we calculate the efficiency rj of the algorithm. 
Qbit 1 represents hereby the "second heat bath" , the 
most elementary bath one can think of. The change of 
the energy expectation value of qbit 1 represents the heat 

Q n = A{E(n)) 1 = Tr{^!p f (n)} - Tv{H 1Pi (n)} (14) 

transferred by the algorithm, with Hi = AEi&z® <to® 0Oj 
Pi(n) and Pf(n) representing the density operator of the 
system before and after the n-th application of the algo- 
rithm without taking the bath contact into account. To 
compute the work W in the system one has to take a 
look at qbit 2 and 3. The energy difference is given by 

A{E(n)) 23 = Tr{H 23 p f (n)} - Tv{H 23Pi {n)} (15) 

with H23 = AE 2 ao® <J Z ® (To + AE3&0® 00® &z ■ If n o work 
was done on the system, the change of the energy ex- 
pectation of this subsystem given by Eq. ()15j) would 
be equal to the heat (fT4"|) with opposite sign, energy 
would only be moved around within the system. Sup- 
posed that I A (2?)23 1 — |A(2?)i| was less than zero, work 
would be extracted from the system, if it was larger 
than zero work would be done on the system. In the 
case that A(E(n)) 23 an d A(E(n))i had different signs, 
\A(E(n)) 23 \ — \ A(E(n))i \ is equivalent to the change of 
the energy expectation value of the entire system, 

W n =A(E(n)) = Tr{Hpi{n)} - Tr{Hpi(n)} 

=Tr{(22 x + H 23 )p { (n)} - Tr^ + H 23 ) Pi (n)} 
= Tr{H 1Pi (n)} -Tr{H 1Pi (n)} 

A(E) 1 

+ Tr{H 2S p ( (n)} - Tr{H 23Pi (n)} . (16) 
s . ' 

A<E) 23 




FIG. 5: The efficiency 771 of the first step n = 1 of the algo- 
rithm 



Now we can define the efficiency r\ n of the algorithm by 
rj n = —Qn/W n - In case of complete relaxation of the 
auxiliary qbits we have thus obtained an analytical result 
for the efficiency r\ n , depicted for the first step for various 
sets of energy splittings of qbit 2 and 3 in Fig. [H In 
the further applications of the algorithm for n — *■ 00 the 
efficiency r\ n goes asymptotically to zero, and, of course, 
always stays below the accounting Carnot efficiency. In 
the region, where the efficiency r\\ is negative, qbit 1 heats 
up instead of cooling down. 

In this article we have presented a cyclic cooling al- 
gorithm. In comparison to the original algorithm intro- 
duced by Boykin et al. in Q we added another quantum 
gate and a bath contact of qbit 2 and 3 using different 
contact times r. The special arrangement of gates and 
bath contacts have made a cyclic application of the al- 
gorithm feasible. For finite bath coupling times r the 
ancilla qbits do not reach their initial temperature and 
also some correlations introduced by the gate operations 
remain within the system. This is not the case for an 
infinite coupling time. Using an adequate r the final 
temperature is lower than the temperature of the orig- 
inal algorithm and approximately reached already after 
less than 10 cycles. 

Furthermore, the cyclic type of algorithm allows for 
a comparison with standard thermodynamical machines 
which always use cyclic processes, e.g., the refrigerator 
cooling a finite reservoir by using mechanical work. Thus, 
we have also computed the used work and the pumped 
heat as well as the efficiency of the algorithm. Finally, 
a comparison to the smallest quantum thermodynamical 
machines (Tol . [l2| comes into reach. However, even those 
non classical algorithms do not break the Carnot limit. 

A further advantage refers to the minimal amount of 
ressources used. The cyclic algorithm always operates 
on the same three qbits. All other improved algorithms 
for lower temperatures are using an exponentially grow- 
ing number of ressources. Of course, in all algorithms 
the final temperature is basically a function of the ini- 
tial one. In order to obtain arbitrary low temperatures a 
combination of different cooling mechanisms is vital. 

We thank M. Henrich, H. Schmidt, H. Schroeder, 
J. Tcifel, G. Rcuther, H. Wcimer, P. Vidal and V. Scarani 
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